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This paper summarises a numerical investigation of phase mixing in time-independent Hamiltonian 
systems that admit a coexistence of regular and chaotic phase space regions, allowing also for low 
amplitude perturbations idealised as periodic driving, friction, and/or white and colored noise. 
The evolution of initially localised ensembles of orbits was probed through lower order moments 
and coarse-grained distribution functions. In the absence of time-dependent perturbations, regular 
ensembles disperse initially as a power law in time and only exhibit a coarse-grained approach 
towards an invariant equilibrium over comparatively long times. Chaotic ensembles generally diverge 
exponentially fast on a time scale related to a typical finite time Lyapunov exponent, but can exhibit 
complex behaviour if they are impacted by the effects of cantori or the Arnold web. Viewed over 
somewhat longer times, chaotic ensembles typical converge exponentially towards an invariant or 
near-invariant equilibrium. This, however, need not correspond to a true equilibrium, which may 
only be approached over very long time scales. Time-dependent perturbations can dramatically 
increase the efficiency of phase mixing, both by accelerating the approach towards a near-equilibrium 
and by facilitating diffusion through cantori or along the Arnold web so as to accelerate the approach 
towards a true equilibrium. The efficacy of such perturbations typically scales logarithmically in 
amplitude, but is comparatively insensitive to most other details, a conclusion which reinforces the 
interpretation that the perturbations act via a resonant coupling. 

PACS number(s): 05.60.+W, 05.40.+j, 51.10.+y, 98.10,+z 



I. INTRODUCTION AND MOTIVATION 

Phase mixing is a fundamental process acting in many- 
particle systems, with important implications for a vari- 
ety of different areas of physics extending from galactic 
dynamics to the physics of charged particle beams 
Pj. In response to external forces and/or self-consistent 
interactions, initially localised collections of particles of- 
ten tend to disperse. Phase mixing is believed to play 
a fundamental role in so-called "violent relaxation," the 
process whereby the stars in an elliptical galaxy evolve 
towards a statistical meta-equilibrium and, under certain 
circumstances, it can play a major role in the dissipation 
of accelerator beams, resulting in substantial - and un- 
desireable - increases in emittance. 

Until recently, most practical considerations of phase 
mixing have assumed, at least implicitly, that the global 
dynamics is integrable or near-integrable, in which case 
initially localised ensembles typically disperse as a power 
law in time. However, there is emerging evidence that, 
in at least some settings relevant to galactic astronomy 
j| and accelerator dynamics 0], one is confronted with 
strongly nonintegrable chaotic dynamics. The manifes- 
tations of phase mixing in chaotic systems are very dif- 
ferent from those in nearly integrable systems @], a fact 
which has led galactic astronomers to speak of a new 
phenomenon called chaotic mixing 0. It is, e.g., obvi- 
ous that localised ensembles of chaotic orbits will typi- 
cally disperse exponentially rather than as a power law 



in time. For this reason, it would seem important to re- 
examine the physics of phase mixing, allowing explicitly 
for systems that admits large measures of chaos. 

Of especial interest are the potential implications of 
chaotic phase space mixing for relaxation towards a sta- 
tistical equilibrium or near-equilibrium. By reformulat- 
ing the dynamics of a time-independent Hamiltonian sys- 
tem as a geodesic flow, e.g., using Maupertuis' Principle 
[Q, some rigorous mathematical results have been de- 
rived. For example, it is known that initially localised en- 
sembles of geodesies on a compact manifold with constant 
negative curvature exhibit a coarse-grained approach, ex- 
ponential in time, towards a uniform population of the 
manifold, i.e., a microcanonical distribution, at a rate A 
that is fixed completely by the magnitude of the curva- 
ture |^] . This implies the existence of a direct connection 
between A and the positive Lyapunov exponent (s) \- If 
the curvature is always negative but not constant, the 
flow becomes significantly more complex, but one can still 
infer an exponential approach towards a microcanonical 
equilibrium at a rate that is bounded from below by the 
least negative value of the curvature |^|. 

However, these results are of limited applicability to 
real physical systems where the curvature can vary in 
sign and, in many cases, may be strictly positive pj. In 
such systems, it would appear that chaos should be 
attributed to parametric instability rather than to nega- 
tive curvature, the division of a constant energy hyper- 
surface into chaotic and regular regions reflecting proper- 
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ties of the local curvature which do, or do not, trigger an 
instability. Within a suitably defined connected chaotic 
phase space region, ensembles may still relax exponen- 
tially towards a microcanonical population p] , but there 
is no reason to expect a one-to-one connection between 
the rate A and the Lyapunov exponents \- 

Additional complications can also arise when the bulk 
dynamics is impacted significantly by phase space struc- 
tures like cantori or an Arnold web |l2|. If a connected 
chaotic phase space region is partitioned by such entropy 
barriers, i.e., phase space obstructions which impede, but 
do not completely block, orbital motions, phase space 
transport can entail a multi-stage process characterised 
by two or more vastly different time scales. In analogy 
with the classical problem of effusion of gas through a 
tiny hole, particles starting on one side of a barrier can 
evolve towards a near-uniform population of the phase 
space region on that side on a time scale t s much shorter 
than the time scale i; on which particles breech the bar- 
rier. On intermediate time scales, t s -C t <C ti, each side 
may be in a near-equilibrium but the densities on the two 
sides will in general be very different. 

Significantly, however, even very small perturbations 
can accelerate phase space transport through such barri- 
ers. It has been long known from the study of maps |l3| ] 
that small random perturbations, i.e., noise, can dramat- 
ically reduce the time scale for diffusion through cantori. 
And similarly, it has been recognised that even very low 
amplitude periodic driving can have important implica- 
tions for phase space transport in systems like plasmas 
[ 14 1 . Recent work on continuous Hamiltonian systems 
[15,[l6| has served to quantify this effect, demonstrating 
that weak noise and/or periodic driving can significantly 
impact phase space transport on time scales much shorter 
than the relaxation time tn on which they significantly 
change the values of quantities like energy, which would 
be conserved absolutely in the absence of perturbations. 
All this would suggest that low amplitude perturbations 
can significantly impact the phase mixing of chaotic or- 
bits, a possibility that has direct physical implications 
for the destabilisation of various structures in a galaxy 
p^t or the defocusing of an accelerator beam Q . 

This paper presents the results of a systematic numer- 
ical investigation of phase mixing for two- and three- 
degrees-of-freedom Hamiltonian systems that admit a 
complex coexistence of regular and chaotic phase space 
regions, allowing also for the effects of low amplitude per- 
turbations idealised as periodic driving and/or friction 
and (in general coloured) noise. 

In principle, one might argue that such an analysis 
is not completely applicable to real many-body systems 
where the dimensionality of the true phase space is much 
larger. However, there are at least three reasons to be- 
lieve that the analysis summarised in this paper should 
find direct physical applications. 

1. In a variety of physical systems, e.g., relatively low 
intensity beams Q or stars in the center of a galaxy or- 
biting close to a supermassive black hole, one expects 



that the bulk properties of the flow will be dominated by 
'external' bulk forces rather than self-interactions; and, 
to the extent that self- interactions are important, one 
might hope to model their effects as friction and noise in 
the context of a Fokker-Planck description. 

2. Recent integrations of gravitationally interacting 
test particles in fixed, time-independent TV-body realisa- 
tions of smooth density distributions have shown that, 
at least for large N, these particles behave in a fash- 
ion closely resembling the motions of particles evolved in 
the corresponding smooth density distribution; and that 
the differences that do exist decrease with increasing N 
p9| Moreover, in terms of both the statistical prop- 
erties of individual orbits and the phase mixing of orbit 
ensembles, discreteness effects are extremely well mim- 
icked by (a suitably defined) Gaussian white noise ]2]J . 

3. Performing honest A^-body integrations of interact- 
ing systems for very large N is not viable using current 
computational resources. However, fully self-consistent 
grid code simulations of intense charged beams have been 
found to exhibit both regular and chaotic phase mixing 
qualitatively similar to what has been found in smooth 
two- and three-dimensional potentials [^p^] . 

White noise appears to be very successful in modelling 
discreteness effects in iV-body systems, where the time 
scale on which the 'random' forces change is extremely 
short. However, real systems are also subjected to other, 
more slowly varying, perturbations which, in many cases, 
one might hope to model as (i) periodic driving and/or 
(ii) coloured noise with a finite autocorrelation time, i.e., 
'random' forces that only change significantly on macro- 
scopic time scales. 

In the context of galactic astronomy, periodic driv- 
ing could reflect the effect of companion galaxies, such 
as the Small and Large Magellanic Clouds orbiting the 
Milky Way. Viewed mathematically, coloured noise is 
a superposition of periodic perturbations combined with 
random phases. One might, therefore, use coloured noise 
to model a high density cluster environment, in which, 
over the course of time, the perturbations acting on any 
given galaxy result from complex interactions with a va- 
riety of different neighbouring galaxies. Alternatively, 
coloured noise might be used to model discrete substruc- 
tures like molecular clouds or globular clusters, where the 
interaction time scale is too long to justify a white noise 
approximation. 

Such perturbations may also influence the dynam- 
ics of charged-particle beams. In a cyclotron or syn- 
chrotron, irregularities in the confining and/or focusing 
magnetic fields would, from the perspective of the beam, 
appear as quasi-periodic perturbations with, perhaps, 
a slowly varying frequency. The conceptual paradigm 
for such beams is the Henon-Heiles potential, V(x, y) = 
\{x 2 + y 2 + 2x 2 y — |y 3 ), although in general detailed 
particle tracking is necessary for determining the effec- 
tive, or 'dynamic', aperture of the machine. The same is 
true conceptually for linear accelerators comprising a pe- 
riodic array of quadrupole and/or solenoid magnets. For 



2 



beams wherein the self-field, i.e., space charge, is active, 
in general the situation is more complicated. Transitions 
in the focusing force can lead to dynamics similar to vi- 
olent relaxation in stellar systems, i.e., the presence of 
a hierarchy of collective modes or, for large transitions, 
correlations analogous to plasma turbulence. Other pro- 
cesses, such as beam generation at the source, or bunch 
compression, may impose irregularities in the form of 
microstructure, i.e., density correlations, on the beam. 
Space-charge forces from these 'clumps' may influence 
individual particle orbits in the manner of colored noise. 

Considerations such as these imply that the micro- 
scopic dynamics of beams and stellar systems are linked, 
a point that, until recently, seems to have been over- 
looked in the literature. The specific mechanism that 
connects these diverse systems is phase mixing. 

Section II of this paper describes the models that were 
considered and the numerical experiments that were per- 
formed. Section III then summarises the results obtained 
for unperturbed Hamiltonian systems. Section IV de- 
scribes how these results are alterred in the presence 
of friction and white noise, which might be expected to 
mimic discreteness effects in a many-body system. Sec- 
tion V considers in turn the effects of periodic driving 
and/or coloured noise with a finite autocorrelation time, 
which might be expected to mimic systematic and/or 
near-random perturbations that cannot be modeled as 
instantaneous kicks. Section VI summarises the princi- 
pal conclusions and comments on potential applications. 



II. THE EXPERIMENTS PERFORMED 

The research described here focused on two- and three- 
degrees-of-freedom Hamiltonian systems of the form 

H= l -{vl+v 2 y +vl)+V{x,y,z). (2.1) 

Three different classes of potential were considered. The 
simplest, which appears to exhibit comparatively generic 
behaviour |2j| , is a straightforward generalisation of the 
two-degrees-of- freedom dihedral potential pifl , allowing 
for three free parameters, a, b, and c: 

\ r { \ 2 2 2 , 1 /„2 , „ 2 , 2\ 2 

V[x, y, z) = -x -y -z + - (x +y + z) 



(ax 2 y 2 + by 1 ?? + cz 2 x 2 ) 



(2.2) 



The second is given as a combination of isotropic and 
anisotropic Plummer potentials [|25(1: 



V(x,y,z) 
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(a 2 + x 2 + y 2 + z 2 ) 1 / 2 



(a 2 + x 2 + b 2 y 2 + c?z 2 ) 1 / 2 ' 



(2.3) 



For generic choices of parameters a, b, c, and m and 'rea- 
sonable' choices of energy E, both these potentials yield 
phase space hypersurfaces admitting a complex coexis- 
tence of both regular and chaotic orbits. 

The third potential is the sum of two integrable con- 
tributions, an anisotropic harmonic potential and an 
isotropic Plummer potential p6[ : 

V(x, y, z) = \{a 2 x 2 + b 2 y 2 + c 2 z 2 ) ~^L=, (2.4) 

with r 2 = x 2 +y 2 +z 2 . This potential is interesting in the 
sense that, provided that M is not too large, chaotic or- 
bits with positive Lyapunov exponent only exhibit strong 
exponential sensitivity for the small fraction of time they 
are at comparatively small r ]2l| . At larger r the motions 
are 'nearly regular.' 

Two-degree-of-freedom systems were obtained trivially 
by considering orbits with initial conditions z = v z = 0. 

Results derived for motion in the unperturbed poten- 
tials were contrasted with motions subjected to two dif- 
ferent classes of perturbations: 

1. Periodic driving. This entailed either allowing for a 
simple sinusoidal driving, resulting in evolution equations 
of the form 

d 2 ri dV . . 

'dt 2= ^dr~ + mUJit (* = x 'y' z )> ( 2 - 5 ) 

or introducing systematic oscillations in one or more of 
the constants in the unperturbed potential, e.g., consid- 
ering a perturbed version of eq. (2.3) with 

a —>■ ao(l + sin tot) or m — + mo(l + sin tot), (2.6) 



V(r,t) (l + Asinut) x V(r), 



(2.7) 



2. Friction and noise. This entailed considering 
Langevin equations of the form 



dt 2 



dV d n 

r\— + F, 



Or 



dt 



(i = x,y,z), (2. 



which incorporate a dynamical friction —rjdri/dt and a 
pseudo-random force Fi , idealised as Gaussian noise with 
zero mean. The friction and noise were typically taken 
to be related by a Fluctuation-Dissipation Theorem 



{F i (t 1 )Fj(t 2 ))=K(t 1 -t 2 )5 ij 
with diffusion constant 



/oo 
dTK(r) = 2? ? e, 
-oo 



(2.9) 



(2.10) 



in terms of a 'temperature' 8 ~ \E\. Three forms of noise 
were considered, namely: 

• white noise, which is delta-correlated in time, with 



3 



K(t) = 2e V S D (r), 



(2.11) 



• coloured noise sampling an Ornstein-Uhlenbeck pro- 
cess, for which the autocorrelation function 



K(r) — ar/Q exp(— a|r|), 



(2.12) 



corresponding to an autocorrelation time t c = 1/a, and 

• coloured noise characterised by a more complex auto- 
correlation function 

K(r) = ^ exp(-a|r|) (l + a\r\ + ^r 2 ) , (2.13) 

for which t c = 2/ a. The white noise autocorrelation 
function (2.10) can be derived from either (2.11) or (2.12) 
in a singular limit a — > oo. White noise was implemented 
numerically using a standard algorithm p8|. Coloured 
noise was implemented using an algorithm developed by 
I. V. Pogorelov as part of his Ph. D. dissertation |p9[ . 

Phase mixing was investigated by tracking the evo- 
lution of localised ensembles of N initial conditions of 
fixed energy E, with 1600 < N < 40000. These were typ- 
ically generated by uniformly sampling a four- or six- 
dimensional phase space hypercube of side AZ, with 
AZ ~ 10 -3 the size of the accessible phase space region. 
The initial conditions were integrated using a standard 
integration scheme which simultaneously tracked the evo- 
lution of a small initial perturbation, periodically renor- 
malised in the usual way |12|, so as to extract estimates 
of the largest (finite time) Lyapunov exponent. 

Two diagnostics were used to probe the evolution of 
each ensemble: 

• lower order moments of the form 
(x A y B z c v x D v^vf) A + B + C + D + E + F < 6, 

(2.14) 

computed in the obvious way; and 

• coarse-grained two-dimensional distribution functions 
like f(x,y,t) or f(x,v x ,t), generated by gridding the or- 
bital data at fixed instants of time into a k x k grid with 
k ranging between 5 and 40. 

The initial divergence of an ensemble was probed 
by tracking the rate at which quantities like the dis- 
persions o~ x or a Va . grow. The approach towards an 
equilibrium or near-equilibrium was probed by tracking 
(i) the convergence of various moments towards time- 
independent 'equilibrium' values and (ii) the evolution 
of coarse-grained distributions towards a nearly time- 
independent form. For example, it is an empirical fact 
that, in each potential, chaotic ensembles typically evolve 
towards a state in which, modulo finite-iV statistics, 
(x) = and o~ x assumes a constant value a x o. The ob- 
vious question, then, is whether the decay of (x(t)} and 
Do~ x (t) = a x (t) — a x o towards zero is power law or expo- 
nential, and, for the case of chaotic ensembles, the degree 
to which the rates associated with this decay correlate 
with the finite time Lyapunov exponents of the orbits. 



Convergence of coarse-grained distributions was han- 
dled analogously Q . Motivated by the visual impression 
that, at least for chaotic ensembles, the orbits usually ap- 
proach a nearly time-independent equilibrium state, the 
last q snapshots of a sequence of p snapshots can be used 
to derive a near-invariant distribution 

1 p 

fniv(Zi, Zj) = - f(Zi,Zj,t q ), (2-15) 

q 

and convergence towards can then be probed by 
computing discretised LP norms 

r,Pf(7 7 +\ — Ifi^j, Zj,t) - fniyjZi, Zj)\ p S \ 

D f ^ z ^>-{ mww ) 

(2.16) 

for p = 1 and 2, with sums extending over the k x k grid. 

III. UNPERTURBED HAMILTONIAN SYSTEMS 

A. Phase mixing of regular ensembles 

Localised ensembles of initial conditions corresponding 
to regular orbits will, when evolved into the future, di- 
verge as a power law in time. In particular, quantities 
like the dispersion a x will typically grow linearly in time 
until they become comparable in magnitude to the size 
of the phase space region sampled by individual orbits in 
the ensemble. 

This power law divergence is easily understood. Mo- 
tion in the x- (or any other) direction is periodic, so that, 
after a period T Xj i — 2h/lu x ^ each orbit i will return to 
its initial value xo,i- In general, however, the periods Tj 
for nearby orbits will be slightly different, so that when 
orbit i has returned to xqa, a slightly displaced orbit j 
with different Tj will not have returned to its initial Xqj . 
It will instead assume a slightly different value Xj which, 
on the average, is slightly further away from xqj than 
the initial Xoj. On the average, the two orbits will have 
moved further apart. (For the special case of a harmonic 
potential, where all the orbits have the same frequencies, 
there is no phase mixing.) 

Let t cr ~ 1/oj denote a typical periodicity, St ~ 1/Su> 
a characteristic 'spread' in periodicities, and R the lin- 
ear 'size' of the configuration space region that is ac- 
cessible to the orbits in the ensemble. The preceding 
argument then suggests that a x should grow in such a 
fashion as to satisfy o~ x /R ~ (Su)t. Presuming, however, 
that Suj ~ u)A/R, where A is the 'size' of the region from 
which the initial conditions were derived, it follows that 

a x (t) - (t/t reg )R, with t reg (A) - {R/A)t cr (3.1) 

or, equivalently, 

a x (t)~ A(t/t cr ). (3.2) 
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That a x grows linearly in time and that, at least 
approximately, a x cx A and t reg cx 1/A are both illus- 
trated in Fig. 1, the first five panels of which track a x (t) 
for ensembles of initial conditions evolved in the two- 
dimensional dihedral potential with a = 1 and energy 
E = 2.0. In each case, a cell of 20000 initial condi- 
tions was constructed, centered about a fixed phase space 
point (xo,yo,v x fl,v y fl). The ensembles were generated 
by setting v y = v y _o, uniformly sampling a configuration 
space square of size A, and, for each choice of x and 
y, solving for v x = v x (E,x,y,v y ) > 0. The five panels 
correspond to values ranging from A = 1.0 x 10 -3 to 
A = 5.0 x 10~ 2 . The sixth panel exhibits on a log-log 
plot the best fit growth time t reg as a function of A. The 
best fit value for the slope is p = —0.95, which implies 
a scaling t reg cx A -0 95 . A dashed line corresponding to 
slope p = —1 is shown for comparison. Given that, for 
this potential and energy, t cr ~ 5 — 10 and R ~ 1, the 
amplitude of t reg (A) is also consistent with eq. (3.1). 

When integrated for much longer times, regular ensem- 
bles typically exhibit a coarse-grained evolution towards 
an invariant equilibrium distribution. This reflects the 
fact that, over very long times, individual regular orbits 
will, in a time-averaged sense, tend to uniformly populate 
the KAM tori to which their motions are restricted. 

This convergence can be quantified by tracking the evo- 
lution of lower order moments for the ensemble, which 
eventually converge towards equilibrium values. Exam- 
ples thereof are exhibited in the first four panels of Fig. 2, 
which were generated from the same ensembles used to 
create Fig. Id and e. These correspond to collections of 
regular 'box' orbits, with the configuration space topol- 
ogy of Lissajous figures. The invariant distributions asso- 
ciated with these ensembles are symmetric with respect 
to reflection about the x- and j/-axes and, as such, have 
moments (x) mv = (y) mv = 0. 

It is evident that (x(t)) is indeed converging towards 
zero, but that the time scale associated with the conver- 
gence is very long. Also evident is the fact that the time 
dependence of (x) is comparatively complex: Strictly 
speaking, the evolution is neither exponential nor power 
law in time; and, even at very late times, t > 4000, one 
can see coherent evolutionary effects reflecting the fact 
that the orbits all have well-defined periodicities. 

One final point is also clear. Just as the size A de- 
termines the time scale on which ensembles disperse ini- 
tially, so also A fixes the time scale on which they evolve 
towards an equilibrium. For the ensemble exhibited in 
panels Figs. 2a and b, A is 2.5 times larger than for 
the ensemble exhibited in Figs. 2c and d, so that, as 
illustrated in Fig. 1, the first ensemble disperses approx- 
imately (2.5) 95 w 2.38 times as fast. However, by re- 
plotting in panels e and f the data for Fig. 2a and b over 
the shorter interval T = (8000/2.38) w 3360, it becomes 
evident visually that the first ensemble also approaches 
an equilibrium approximately 2.38 times as fast. 



B. Phase mixing of chaotic ensembles 

Initially localised ensembles of chaotic orbits evolve 
very differently. Considerable variability can be observed 
but, overall, chaotic ensembles tend to disperse exponen- 
tially at a rate A = t ~ ha that is comparable to the largest 
positive Lyapunov exponent x- Thus, e.g., 

a x ~ A exp(*/t cha ) with t cha ~ x~\ (3-3) 

with A the 'size' of the cell of initial conditions. This 
exponential growth implies that chaotic phase mixing is 
typically much faster than regular phase mixing. 

The physics of chaotic phase mixing is very different 
from the phase mixing of regular ensembles. Aperiodic 
chaotic orbits exhibit exponentially sensitive dependence, 
so that two nearby initial conditions separated by some 
phase space distance 5Zq will, when evolved into the fu- 
ture, typically diverge exponentially at a rate that is set 
by a Lyapunov exponent: SZ(t) ~ SZq exp(xt). Apply- 
ing this principle to the ensemble as a whole leads imme- 
diately to eq. (3.3). 

The appropriate x entering into eq. (3.3) is a, finite time 
Lyapunov exponent [ |30"[ appropriate for orbits in the en- 
semble for the interval associated with the exponential 
divergence, not the asymptotic Lyapunov exponent, as 
defined in a t — > oo limit, which can be very different. 
Two different ensembles that rapidly spread to occupy 
the same phase space region might be expected to have 
comparable x's and to mix at comparable rates. How- 
ever, it is possible for different ensembles, each in prin- 
ciple able to access the same phase space region, to have 
different finite time x's and to mix at very different rates. 

If the dispersal of the ensemble be probed in a generic 
direction, t'Z ia will typically be comparable to the largest 
positive (finite time) Lyapunov exponent. In general, a 
conservative nonintegrable potential in D dimensions has 
D — 1 unequal positive Lyapunov exponents, correspond- 
ing to D — 1 different unstable directions. However, the 
rate at which an ensemble disperses in a specified phase 
space direction will be dominated by the Lyapunov ex- 
ponent associated with the most unstable direction with 
a nonzero projection in that specified direction; and that 
most unstable direction will typically have a nonzero pro- 
jection in a generic, 'randomly chosen' direction pjj . 

It is, however, possible for certain directions to mix 
much more rapidly than others. As an extreme exam- 
ple, consider a three-dimensional potential where motion 
in (say) the z-direction is integrable and decoupled from 
chaotic motion in the two orthogonal directions. In this 
case, quantities like o~ z and a Vz will exhibit behaviour 
appropriate for regular phase mixing whereas quantities 
like a x or a Vy will in general exhibit the exponental be- 
haviour predicted by eq. (3.3). 

Individual chaotic orbits evolved in a generic poten- 
tial can be impacted significantly by such topological 
obstructions as cantori and/or Arnold webs ]12| , which 
serve as entropy barriers to impede - but not block - 
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motion between different portions of the same connected 
phase space region. If, however, the motions of individual 
members of an ensemble are impeded, the phase mixing 
of the ensemble as a whole will in general be alterred. 

As a simple example, consider an ensemble of initial 
conditions corresponding to chaotic orbits which, at early 
times, are trapped near a regular island by cantori. The 
largest finite time Lyapunov exponents for such 'sticky' 
p2[ chaotic orbits will in general be substantially smaller 
than the exponents for 'wildly chaotic' orbits not trapped 
near regular islands ]33[ ] ; and this implies that the ensem- 
ble will disperse more slowly than a 'wildly chaotic' en- 
semble. Even more significant, however, is the fact that, 
until the orbits have escaped through the cantori, quanti- 
ties like a x may remain much smaller in magnitude than 
the asymptotic value a Xl00 towards which the ensemble 
converges at late times. 

Representative data for a 'wildly chaotic' ensemble are 
exhibited in the left panels of Fig. 3, the upper curves in 
which exhibit a x , a y , and a z for an ensemble evolved in 
the potential (2.4) with a = 1.95, b = 1.50, c = 1.05, and 
M = 0.0316. These curves show considerable short time 
structure but it is evident that, overall, the dispersions 
grow exponentially until they saturate on a scale com- 
parable to the size of the accessible configuration space 
region. The exponential growth becomes especially ap- 
parent if the data are smoothed via box-car averaging. 
The lower curves in these panels, translated downwards 
by a distance lncr = 3.0, illustrate the effect of averaging 
the raw data, recorded at intervals St = 0.25, over 20 
adjacent points. It is apparent visually that the growth 
rates associated with <j x and a y are essentially identical. 
The behaviour in the z-direction is noticeably different, 
but the overall growth rate remains comparable. 

Fig. 4 exhibits dispersions for chaotic ensembles 
evolved in the potential (2.2) with a = b = c = 1. 
In each case, data were recorded at intervals St = 0.5, 
and the curves generated by box-car averaging over 5 
adjacent points. Panels a and b present, respectively, 
smoothed configuration and velocity space dispersions 
for a 'wildly chaotic' ensemble with E = 1.0, with the 
y- and z-components translated downwards by, respec- 
tively, lncr = 5 and 10. The dots surrounding the top 
curves represent unsmoothed data. It is evident that, for 
this ensemble, all six dispersions grow exponentially at 
comparable rates. Panel c presents configuration space 
dispersions for a 'sticky' ensemble with E — 4.0. In this 
case, one still identifies an initial exponential growth, but 
that growth slows for lncr > —1.0 or so, and the even- 
tual approach towards a near-invariant cr is accompanied 
by oscillations not observed in most 'wildly chaotic' en- 
sembles. Panel d exhibits exhibits configuration space 
dispersions for a particularly unusual example, in which 
dispersal in the z-direction is much slower than in the x- 
and y-directions. This anomalous behaviour reflects the 
fact that the cell of initial conditions was concentrated 
near the x — y plane, where dV/dz = 0. 

Atypical behaviour is especially common for potentials 



like (2.4), where chaotic orbits behave in a nearly regular 
fashion 'most of the time.' It is, e.g., easy to identify en- 
sembles which, at least initially, avoid the 'highly chaotic' 
central region and, as such only mix on comparatively 
long time scales. And, similarly, it is easy to construct 
ensembles which, because of the (initial) orbital orienta- 
tions in the central regions, are considerably more chaotic 
in some phase space directions than in others. 

As a clear example, one can contrast the data for the 
'wildly chaotic' ensemble exhibited in the left panels of 
Fig. 3 with the right hand panel, which exhibits data 
for a different, much less chaotic, ensemble evolved in 
the same potential with the same energy. For the left 
ensemble, the mean value of the largest finite time Lya- 
punov exponent at t = 256 is (\) ~ 0.055; for the right 
ensemble, (x) ~ 0.022. It is clear that, overall, the dis- 
persions grow more rapidly for the wildly chaotic ensem- 
ble. Overlaying plots for the two ensembles also reveals 
another significant point: Although the sticky ensemble 
originally diverges more slowly, the rates of divergence 
become comparable when hxa x > —2, the growth rate 
for the sticky ensemble having increased substantially. 
This increase reflects the fact that, by this time, most of 
the 'sticky' orbits have become unstuck. 

The initial exponential divergence eventually saturates 
and, if the ensembles are probed for longer times, one 
typically observes a coarse-grained evolution, again expo- 
nential in time, towards an invariant, or near-invariant, 
distribution fi n iv> a statistical (near-)equilibrium. In 
particular, reduced distributions f(Zi,Zj) and various 
moments typically converge exponentially towards (near- 
equilibrium values, so that, e.g., 

Da x (t) = a x (t) - a xfi ~ exp(-Ai). (3.4) 

This near-invariant fx n i V appears to correspond to a uni- 
form population of those portions of the phase space 
which are easily accessible to the orbits, i.e., blocked by 
neither conservation laws nor significant entropy barriers. 
If, e.g., one considers an ensemble of orbits, all with the 
same energy Eq, evolved on a phase space hypersurfacc 
that is almost completely chaotic and where topological 
obstructions are unimportant, the distribution towards 
which the ensemble evolves is well approximated by a 
microcanonical distribution fi oc Sd{E — Eq). |3~f| 

Just as, in many cases, the ensemble initially dis- 
perses at roughly the same rate in different directions, 
so also it is often true that the approach towards a near- 
equilibrium proceeds at comparable rates in different di- 
rections. Indeed, at least for 'wildly chaotic' ensembles, 
the rate A associated with the approach toward a near- 
equilibrium appears typically to be comparatively insen- 
sitive to the diagnostics which one chooses to probe. 

As probed by lower order moments, the rate of conver- 
gence is insensitive to choice of phase space direction, so 
that, e.g., Dhi(x), Dla(y), and Dhi(v z ) typically con- 
verge towards zero at comparable rates. The order of the 
moment is also unimportant, so that, e.g., the conver- 
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gence rate associated with D In (x p ) is relatively insensi- 
tive to the choice of p, at least for p < 6. Cross moments 
also evolve similarly, so that, e.g., the rates associated 
with Dln(xy} and Dln(xv x ) are comparable. Analo- 
gously, the convergence rates associated with different 
reduced distributions like D p f(x,y) and D p f(x,v x ) are 
comparable to one another and relatively insensitive to 
the level of coarse-graining, i.e., the choice of k for the 
k x k binning, at least for 5 < k < 40. And, perhaps most 
strikingly, the rates associated with D p f are comparable 
to the rates associated with lower order moments. 

Figs. 5 and 6 exhibit the evolution of a variety of dif- 
ferent moments and reduced distributions for two differ- 
ent ensembles evolved in the potential (2.2) with a = 
b = c = 1. Fig. 5 was generated from a 'wildly chaotic' 
ensemble with E = 2.0. Fig. 6 was generated from a 
'sticky' ensemble with E = 4.0. It is evident that all the 
curves in Fig. 5 exhibit an initial exponential decrease; 
and it is also apparent that the slope associated with dif- 
ferent moments and different reduced distributions are 
comparable in magnitude. (The approach towards non- 
zero nearly constant values at later times is a reflection 
of finite N statistics.) By contrast, many of the curves 
in Fig. 6 exhibit significant deviations from the linear 
behaviour associated with an exponential decrease; and, 
even when the data are well fit by an exponential, the 
rates can differ appreciably. 

The extent to which the moments for these ensem- 
bles are, or are not, comparable can be gauged visually 
from the first two panels of Fig. 7, which record best 
fit rates A for 32 representative quantities. The first 12 
quantities, represented by diamonds, correspond to odd- 
powered moments (x), (y), (z), (v x ), (v y ), (v z ), (x 5 ), 
(y 5 ), (z 5 ), (v x ), ("Vy), and (v z ) . The next six, repre- 
sented by squares, correspond to even-powered moments 
(x 2 ), (y 2 ), (z 2 ), (v 2 ), (v 2 ), and {v 2 z ). Following these are 
mixed moments (xy), (yz), and (zx), represented by x's, 
and (xv x ), (yv y ), and (zv z ), represented by +'s. The 
last eight points, represented by triangles, correspond to 
the quantities D 1 f(v x ,v y ), D 1 f(x,v x ), and D 1 f(x,y), 
computed for k = 40 with an L 1 norm, D 1 f(x,y) for 
k = 30,20, and 10, and the quantities D 2 f(x,y) and 
D 2 f(v x ,v y ) computed with an L 2 norm. 

For the 'wildly chaotic' ensemble used to generate 
Fig. 7a, all the moments and most of the reduced dis- 
tributions converge at similar rates, comparable to the 
mean rate associated with all the moments, which is rep- 
resented by the horizontal line. The largest deviations 
are observed for the L 2 distances D 2 f, where the conver- 
gence rates are appreciably larger. The 'sticky' ensemble 
in Fig. 7b is more complex, with the best fit rates exhibit- 
ing a much larger variability. There are, however, certain 
regularities. For example, moments like (xy), (yz), and 
(zx), which probe two different configuration space di- 
rections, converge substantially more slowly than (xv x ), 
(yv v ), and (zv z ). Overall, it would appear that the mo- 
ments divide empirically into two classes, corresponding 
to faster and slower convergence rates. The two vertical 



lines represent the mean rates for the moments with rates 
greater than, and less than, A = 0.042. 

The rates A associated with the approach towards equi- 
librium are invariable smaller than the typical value of 
the largest (finite time) Lyapunov exponent x- However, 
one often observes strong correlations between the values 
of A and x, larger \ correlating with a more rapid ap- 
proach towards a near-equilibrium. Significantly, though, 
this correlation does not appear universal. 

An example of the connection between A and \ is ex- 
hibited in panel (c) of Fig. 7, which exhibits the ratio 
X/x as a function of energy E for some typical 'wildly 
chaotic' ensembles evolved in the potential (2.2). Dia- 
monds represent values generated by averaging over the 
rates associated with (x) , (y) , and (z) , whereas triangles 
represent values generated from Df(x,y), Df(y,z), and 
Df(z,x). In each case, the error bars reflect variances 
associated with the mean values. The ratio X/x seems 
to assume a nearly constant value -~ 0.2, independent 
of energy. However, this ratio exhibits significantly more 
variability for 'sticky chaotic' ensembles, although it does 
appear to remain smaller than unity. 

The preceding examples illustrate the fact that, for 
'wildly chaotic' ensembles, one can speak meaningfully 
of a coarse-grained approach, exponential in time, to- 
wards a near-invariant distribution on a time scale t n i V 
that is comparable to, but somewhat longer than, the 
characteristic Lyapunov time X • For 'sticky' ensem- 
bles the evolution will be somewhat slower, and sizeable 
deviations from a purely exponential evolution may be 
apparent, but again there is a roughly exponential ap- 
proach towards a near-invariant distribution. 

In neither case, however, is there a guarantee that this 
near-invariant distribution corresponds to a true equi- 
librium. Because of entropy barriers like cantori or an 
Arnold web, orbits may be impeded temporarily from ac- 
cessing various phase space region which are in principle 
accessible, so that the ensemble may evolve initially to- 
wards a distribution that involves a nearly uniform pop- 
ulation of only part of the accessible phase space hyper- 
surface. Only over considerably longer time scales is a 
true equilibrium achieved. 

This is illustrated in the first two panels of Fig. 8, 
which were generated from a collection of 20000 initial 
conditions with energy E = 6.0 evolved in the two- 
dimensional dihedral potential with a = 1. The members 
of the ensemble, corresponding initially to wildly chaotic 
orbits, rapidly approached a near-invariant distribution, 
the form of which can be inferred from Fig. 8a, which ex- 
hibits the spatial coordinates of all the orbits at t = 36.0. 
This distribution is similar to, but clearly distinct from, 
the apparent true invariant distribution which is only ap- 
proached on a significantly longer time scale. This true 
invariant distribution, illustrated in Fig. 8b by the spa- 
tial coordinates at t = 720, is clearly more nearly uniform 
than the earlier near-invariant distribution. In particu- 
lar, the distinctive 'x-shaped' region of higher concentra- 
tion has been blurred as orbits diffused away from the 
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diagonals to occupy regions that were avoided originally. 

IV. THE EFFECTS OF FRICTION AND WHITE 
NOISE 

For sufficiently large amplitudes and/or over suffi- 
ciently long time scales, friction and white noise will 
dramatically effect the evolution of orbit ensembles since 
energy is no longer conserved. For 'typical', e.g., addi- 
tive, white noise, the energies of individual orbits only 
change significantly on a relaxation time scale tn = T) , 
as the ensemble evolves towards a canonical distribution, 
/ oc exp(— E/Q), although various sorts of multiplicative 
noise can dramatically reduce £r [p5[ . The objective here 
is to assess the effects of friction and noise over time scales 
sufficiently short that the energies of individual orbits are 
almost conserved. Interest, therefore, will be restricted 
largely to effects proceeding on time scales <C rf l ■ 

Because random perturbations typically push nearby 
trajectories apart, friction and noise tend generically to 
increase the efficacy of phase mixing for both regular and 
chaotic ensembles. 

Phase mixing of unperturbed regular orbits arises be- 
cause nearby orbits typically oscillate with unequal fre- 
quencies. Perturbing the orbits by friction and white 
noise gives rise to an additional phase mixing that is very 
different in origin. Friction and noise act to induce a ran- 
dom walk in phase space, so that an unperturbed orbit 
and a noisy orbit with the same initial condition, or two 
different noisy realisations of the same initial condition, 
will diverge, at least initially, as t 1 / 2 . Consider, e.g., a 
harmonic potential, where one can proceed analytically 
p8[ . Assuming that the natural frequency lu ^ ?/, one 
can average over oscillatory behaviour proceeding on a 
time scale uj^ 1 to conclude that, for t <C 77 -1 , an ensem- 
ble generated from different noisy realisations of the same 
initial condition will disperse in such a fashion that, e.g., 

uj 2 Sx 2 = 2G??t. (4.1) 

This is a comparatively weak effect, Sx growing on the 
same time scale as the spread in energies which satisfies 
5E 2 /E = 2Qrjt, with E the initial energy. 

For more generic potentials, where unperturbed orbits 
oscillate with unequal frequencies, the evolution can be 
more complex. Friction and noise can be viewed as con- 
tinually deflecting orbits from one periodic trajectory to 
another, but orbits with different periodicities are sus- 
ceptible to the linear phase mixing described in the pre- 
ceding section. One point, however, is clear. For the case 
of regular orbits, friction and noise can have significant 
effect on the bulk evolution of an ensemble only on the 
long time scale 77" 1 associated with changes in energy. 

Nevertheless, friction and noise can have noticeable ef- 
fects over shorter time scales by 'fuzzing out' shorter scale 
structure. Because regular orbits are periodic, a plot 
of a x for an unperturbed regular ensemble will typically 



exhibit considerable structure superimposed on an aver- 
age linear growth law. The introduction of friction and 
noise destroys this exact periodicity and, consequently, 
can help suppress this structure. At a microscopic level, 
this 'fuzzing' is associated with the fact that, with the 
introduction of noise, Liouville's Theorem is no longer 
applicable, so that it becomes possible for orbits to self- 
intersect. Friction and noise can also accelerate the evolu- 
tion towards a 'well-mixed' state in which the basic sym- 
metries of the potential are manifest, so that, e.g., for the 
potentials described in Section II, (x) = (y) = (z) = 0. 
However, this is comparatively slow effect, which only 
proceeds on the time scale tn = r\~ x associated with sig- 
nificant changes in the energies of individual orbits. 

As an example, consider Fig. 9, which exhibits In |(a;)| 
for an ensemble of regular orbits with initial energy 
E = 2.0 evolved in the potential (2.2). The first panel, 
generated from unperturbed orbits, shows no evidence of 
evolution towards (x) = 0. The remaining three panels, 
generated from orbits evolved with 9 = E and ?] = 10~ 7 , 
1CP 6 , and 10~ 5 , show clear evolutionary effects, but only 
on a timescale comparable to t^. In each case, the energy 
dispersion is well fit by a linear growth law, 

a% = A® V E t, (4.2) 

with Eq the initial energy and A a constant of order unity, 
so that, e.g., at time t — 512, the energy dispersions for 
the three panels were, respectively, oe = 0.032, 0.101, 
and 0.318. It is thus evident that, by the time has 
begun to evidence a significant systematic decrease, the 
orbital energies have manifested appreciable changes, so 
that one can no longer envision an ensemble of orbits 
restricted, even approximately, to a fixed constant-energy 
hypersurface. 

Noise acts to enhance chaotic phase mixing in a very 
different way. Because of their exponential sensitivity, 
perturbed and unperturbed orbits with the same initial 
conditions, or two different noisy realisations of the same 
initial condition, will typically diverge at a rate set by an 
appropriate finite time Lyapunov exponent, but with a 
prefactor that depends on the amplitude of the pertur- 
bation. In particular, one has J37J 

8r~J-±exp(t/tcka), (4-3) 
V uj 6 

with t c ha ~ X an( l T ~ 1/w a characteristic orbital time 
scale. Noise acts as a 'seed' to push apart two orbits with 
the same initial condition but, once separated, their sepa- 
ration will grow exponentially. The form of the prefactor 
is easily understood: Noise acts as a 'random' process, 
so that its efficacy scales as the square root of the ampli- 
tude of the perturbation, i.e., oc (O??) 1 / 2 . The remaining 
dependence on u> follows from dimensional analysis. 

The implication for chaotic mixing is obvious: If the 
size A of the initial ensemble is extremely small or if 
the amplitude of the noise is comparatively large, noise 
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can have a substantial effect on the rate at which quan- 
tities like o x grow, even for 'wildly chaotic' ensembles. 
This could, for example, be of practical importance for 
charged particle beams, where the aim is to minimise the 
growth of omittance in a (hopefully) very narrow beam. 
If, however, A is not that small, noise will have a com- 
paratively minimal effect on the rate at which 'wildly 
chaotic' ensembles mix, although it can prove important 
by helping to smooth out short scale structure |56| . 

For 'sticky' ensembles, weak noise can again smooth 
out short scale structure. Even more importantly, how- 
ever, it can make the ensemble less 'sticky'. Indeed, if the 
noise is sufficiently strong, differences between 'sticky' 
and 'wildly chaotic' behaviour can be completely erased, 
an otherwise 'sticky' ensemble evolving in the fashion 
that appears 'wildly chaotic' 

Both these effects are illustrated in Fig. 10, which ex- 
hibits <r x for two different chaotic ensembles evolved in 
the dihedral potential (2.2) in the presence of noise of 
variable amplitude ranging from 77 = to 77 = 10~ 4 . In- 
spection of the left hand panels, generated from a 'wildly 
chaotic' ensemble with E = 2.0, shows that the noise has 
only a comparatively minimal effect on the growth of a x . 
By contrast, for the 'sticky' ensemble with E = 4.0 ex- 
hibited in the right hand panels, noise with amplitude as 
small as 77 = 10~ 5 is sufficient to largely obliterate the ir- 
regular structure that, in the absence of noise, is evident 
at very early times. For noise with amplitude as large 
as ?7~ 4 , all structure has completely disappeared. It is 
evident visually that the top two curves in the Figure, 
corresponding to 77 — 10~ 4 , both look 'wildly chaotic'. 

Weak noise can prove even more important for the 
convergence of an ensemble towards an invariant or near- 
invariant distribution. Even for wildly chaotic ensembles, 
noise acts to 'fuzz out' orbits on short scales, so that 
the rate of approach towards a near-invariant distribu- 
tion can be accelerated significantly and, for 'sticky' en- 
sembles, this effect can be especially dramatic, the noise 
serving also to suppress structure that is otherwise man- 
ifested in quantities like ln|(x)| or Df(x,y). Examples 
thereof are illustrated in Fig. 11, which exhibits Df{x, y) 
for the same two ensembles used to generate Fig. 10. 

It is clear from even a cursory examination of Figs. 10 
and 11 that the effects of noise on the initial disper- 
sal of an ensemble and its eventual approach towards a 
near-invariant distribution must exhibit at most a weak, 
roughly logarithmic dependence on 77. This is confirmed 
by Fig. 12, which plots the rate A with which f(x,y) 
approaches a near-invariant distribution as a function of 
logic) 7 /" 1 - For V < 10~ 6 - 5 or so, i.e., t R > 10 65 , the 
noise has a comparatively minimal effect. For larger 77 
and smaller 4r, however, there is an obvious logarithmic 
dependence. The left-most two points in the Figure, gen- 
erated for the shortest tR, correspond to slightly higher 
values of A than would have been expected from a simple 
extrapolation from the remaining points. This, however, 
easily understood, reflecting as it does the fact that, for 
such very large values of 77, energy nonconservation is 



important even on time scales t <C 50. 

This logarithmic dependence is reminiscent of the log- 
arithmic dependence on r\ observed for noise-enhanced 
diffusion through cantori [lq| . Consider, e.g., multiple 
noisy realisations a single initial condition corresponding 
to an orbit trapped near a regular island by cantori. In 
this case, one finds typically that, once the orbits have 
dispersed to fill the easily accessible regions not blocked 
by the cantori, they will begin to leak through the can- 
tori in a fashion that is well modeled by a Poisson process 
with a rate A that scales logarithmically with amplitude. 

Finally, it should be observed that, by facilitating dif- 
fusion through entropy barriers, even weak noise can dra- 
matically accelerate the approach towards a true invari- 
ant distribution. This effect is illustrated in the bottom 
two panels of Fig. 8, which exhibit orbital data for the 
same initial conditions as the top panels but now evolved 
in the presence of noise with = E and 77 = 10~ 6 and 
10~ 5 . The data in both these panels were extracted from 
noisy orbits at time t = 36.0. However it is evident vi- 
sually from a comparison with Fig. 8a and b that they 
more closely resemble unperturbed orbits at the much 
later time t = 720 than at t — 36.0. 



V. THE EFFECTS PERIODIC DRIVING AND 
COLOURED NOISE 

Provided that the autocorrelation time t c is short com- 
pared with the natural time scale t cr for the orbits in 
some ensemble, its precise value is unimportant, and 
coloured noise has nearly the same effect, both quali- 
tatively and quantitatively, as white noise with the same 
diffusion constant D. However, when t c becomes com- 
parable to t cr , the effects of the noise grow substantially 
weaker; and, for t c 3> t cr , the noise has a comparatively 
minimal effect. Overall, for t c ^> t cr the effects of the 
noise exhibit a roughly logarithmic dependence on t c . 

Fig. 13 demonstrates how the evolution of the un- 
perturbed ensemble exhibited in Figs, lib and 12b is 
changed by incorporating Ornstein-Uhlenbeck noise with 
= 4.0, 77 = 10~ 5 , and different autocorrelation times 
ranging between t c = 0.01 and t c = 10.0. For t c < 0.33 
the precise value of t c matters little, the evolution of Df 
being virtually identical to what was observed for the 
case of white noise with the same amplitude. Alterna- 
tively, for t c > 3.33, the effect of the noise is drastically 
reduced, a x and Df behaving very nearly as if no noise 
were present. Fig. 14 exhibits analogues of Fig. 13, gen- 
erated now for orbits evolved in the presence of the more 
complicated noise with autocorrelation function given by 
eq. (2.13). It is clear visually that the two noises have 
virtually identical effects. 

Figs. 15 a and b quantify this accelerated evolution 
by tracking the rate A associated with the approach of 
f{x,y) towards a near-invariant distribution. The data 
points were obtained for a wildly chaotic ensemble with 







initial energy E = 2.0 again evolved in the dihedral po- 
tential, but now allowing for coloured noise with = 1.0, 
7/ = 10~ 5 , and variable t c . The top panel was gen- 
erated for orbits evolved in the presence of Ornstein- 
Uhlcnbeck noise; the lower for orbits evolved with noise 
characterised by the autocorrelation function (2.13). For 
( c « 1, the precise value of t c is largely irrelevant and 
A assumes, at least approximately, the value that ob- 
tains for white noise. Alternatively, for t c ^> 1 the rate 
A assumes a near-constant value close to that associated 
with the evolution of an unperturbed ensemble. For in- 
termediate values, A is a decreasing function of the auto- 
correlation time which scales logarithmically in t c . Also 
evident is the fact that the two differences noises have 
virtually identical effects. That the scaling in t c is again 
logarithmic and that the precise form of the noise seems 
unimportant is again reminiscent of what has been ob- 
served for noise-enhanced diffusion through cantori [ p"5| . 

These effects are easily understood qualitatively, given 
the expectation that noise impacts individual orbits via 
a resonant coupling |jl5). Most of the power in an orbit, 
either regular or chaotic, is concentrated at frequencies 
to r-j 2tt /t cr , with t cr the natural orbital time scale, and an 
efficient coupling requires noise with substantial power at 
comparable frequencies. Delta-correlated white noise has 
a flat power spectrum and, as such, can couple to more 
or less anything. However, the introduction of a finite 
autocorrelation time suppresses high frequency power. If, 
in particular, t c 3> t cr there is essentially no power at the 
natural frequencies associated with the orbits, so that the 
noise has a comparatively minimal effect. In this case, 
the only coupling arises via higher order harmonics. 

This interpretation assumes implicitly that the effects 
of noise are not associated with a few especially large 
kicks, but result rather from a continuous wiggling trig- 
gered by 'typical' kicks. This seems reasonable when con- 
sidering noise as a source to 'fuzz' out trajectories within 
a single, easily accessible phase space region; and an anal- 
ysis of diffusion through cantori and/or along an Arnold 
web jy| suggests that, at least for these sorts of entropy 
barriers, this is also the case. By continually wiggling the 
orbits, noise helps them 'hunt' for holes in cantori. 

This suggests further that the form of the noise should 
be largely irrelevant in determining its overall effect on 
phase space transport. At least within the class of Gaus- 
sian noise, this appears to be true. It may be possible to 
fine-tune the noise to have especially large or small effects 
but, generically, the details appear to be unimportant. If, 
e.g., the additive white noise described in the preceding 
section is replaced by multiplicative noise by making rj 
a function of the phase space coordinates, very little is 
changed. Thus, if the constant 77 — rj considered there 
is replaced by a new ij — rjo(v/ (v)) ±2 , with (v) the aver- 
age speed of the unperturbed orbits, one finds that the 
rates at which moments and coarse-grained distributions 
evolve towards a near-invariant distribution remain es- 
sentially unchanged. Similarly, the form of the colour 
is comparatively unimportant. For example, noises with 



the autocorrelation functions given by eqs. (2.11) and 
(2.12) yield very similar effects. The presence or ab- 
sence of an accompanying dynamical friction also appears 
largely irrelevant. 

That noise acts on orbits via a resonant coupling can be 
corroborated by studying how, when subjected to time- 
dependent perturbations, the unperturbed energies of the 
unperturbed orbits are changed: Even though the effects 
described here are not associated with changes in energy 
per se, the degree to which energy is changed provides 
a simple probe of the strength of the coupling. Since 
coloured noise can be interpreted as a superposition of 
periodic drivings with different frequencies, it is useful to 
decompose the noise into its basic constituents by study- 
ing individually the effects of periodic driving with dif- 
ferent frequencies to. The results of such an investigation 
is exhibited in the top two panels of Fig. 16. The data in 
Fig. 16a, produced for the 'sticky' ensemble used to gen- 
erate Figs. 10 and 11, was generated by reintegrating the 
orbits in the presence of a perturbation of the form given 
by eq. (2.5), with amplitude A x = A y = A z = A = 10~ 2 
but variable lo x = uj y = io z = u>. Each integration pro- 
ceeded for a total time t = 256, and the energy dispersion 
cje was recorded at intervals St = 0.25. The quantity 
plotted is the maximum value assumed by oe during the 
integration at times t > 10.0 (after the decay of any tran- 
sients). Fig. 16b shows analogous data, generated for the 
same ensemble of initial conditions but now integrated in 
a perturbed potential of the form given by eq. (2.7). 

It is evident in both cases that, for low frequencies, oe 
assumes a near-constant value, that there is a sharp in- 
crease for intermediate frequencies, and that, for higher 
frequencies, <je is again comparatively small. For per- 
turbations satisfying eq. (2,7), oe oc 1/lo for large lo; for 
perturbations satisfying eq. (2.5), oe is nearly indepen- 
dent of lo for high frequencies. For perturbations sat- 
isfying eq. (2.7), the detailed response at intermediate 
frequency exhibits a comparatively complex dependence 
on u; for those satisfying eq. (2.5), the dependence on lo 
is smoother. For high and low lo, the time dependence of 
<TE(t) can exhibit considerable variability; but, for inter- 
mediate frequencies where the response to the perturba- 
tion is strong, one observes invariably a response well fit 
by a growth law a E oc t, the same growth law associated 
with (both white and coloured) noise. 

The connection with properties of the unperturbed or- 
bits is apparent from Fig. 16c, which exhibits Fourier 
spectra generated from an evolution of the unperturbed 
ensemble for a time t = 256 with data recorded at inter- 
vals St = 0.05. The quantities |ar(o;)| (solid curve) and 
|j/(w)| (dashed) plotted there were generated by comput- 
ing the quantities |a;i(a;)| 2 and |j/i(w)| 2 for each orbit in- 
dividually, and then constructing, e.g., a composite 

N 

IsMf^kiMI 2 . (5.1) 

i=l 

It is obvious that these power spectra share certain fea- 
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tures with the plots of ue- One again sees a flat depen- 
dence on w at low frequencies, substantial power at in- 
termediate frequencies, and a power law decrease, oc l/u>, 
for sufficiently large frequencies. The peak frequencies in 
Fig. 16a coincide with the peak frequencies associated 
with the Fourier spectra, which suggests a simple reso- 
nant coupling. The peak frequencies in Fig. 16b occur 
at somewhat higher frequencies. This suggests couplings 
through harmonics may be more important in this case, 
which is hardly surprising, given the more complex time- 
dependence entering into (2.7). 

Given these results, it is easy to predict the effects 
of periodic driving on the evolution of orbit ensembles. 
Since coloured noise is a superposition of drivings with 
different periodicities, one would expect the same loga- 
rithmic dependence on amplitude; and, by analogy with 
the case of noise, one would only expect a significant re- 
sponse if the frequency u) associated with the driving is 
comparable to the natural frequencies of the orbits that 
are being perturbed. However, one might anticipate that 
quantities like the rate A associated with the approach to- 
wards a near-invariant distribution would exhibit a more 
irregular dependence on u) than the relatively smooth 
variation of A(t c ) associated with coloured noise. All 
these predictions are in fact correct. 

Fig. 17 contrasts the evolution of an unperturbed en- 
semble with E = 4.0 integrated in the two-dimensional 
dihedral potential with the same ensemble integrated in 
the presence of periodic driving of the form (2.7) with 
amplitude A = 10 -15 and variable frequencies ranging 
between oj = 0.1 and oj = 30.0. The left panel exhibits 
the dispersion a x both for raw data recorded at intervals 
St = 0.05 (dots) and data smoothed via a box-car aver- 
age over 20 adjacent points (solid lines). The right panel 
exhibits Df(x,y). Two points are apparent. First, it is 
clear that although the driving has a significant effect for 
driving frequencies oj = 3.0 and oj = 10.0, the higher and 
lower frequencies have a comparatively minimal effect. 
This is completely consistent with the fact that, as is ev- 
ident from Fig. 15 (b), the strongest coupling between 
the driving and the unperturbed orbits is for energies 
between roughly 2.0 and 20.0. The second point is that, 
for oj — 3.0 and 10.0, the driving impacts a x both by 
increasing the overall rate of increase and by suppressing 
the prominent wiggles which are evident for the unper- 
turbed ensemble and for the ensembles perturbed with 
higher and lower frequencies. 

VI. DISCUSSION 

This paper has summarised an investigation of phase 
mixing in 'realistic' Hamiltonian systems which admit 
both regular and chaotic orbits, focusing also on the ef- 
fects of low amplitude perturbations modeled as periodic 
driving and/or friction and (in general coloured) noise. 

Localised ensembles of regular initial conditions will, 



when integrated into the future, diverge linearly in time 
so that, e.g., the dispersion a x ~ A(t/t cr ), with A the 
'size' of the region sampled by the initial conditions and 
t cr a characteristic orbital time scale. By contrast, ini- 
tially localised ensembles of 'wildly' chaotic orbits diverge 
exponentially at a rate which is typically comparable to 
the magnitude of the largest finite time Lyapunov expo- 
nent for orbits in the ensemble. 'Sticky' chaotic orbits 
can exhibit more complex behaviour, but again tend to 
diverge exponentially. The prcfactor associated with the 
exponential growth of quantities like a x again exhibits a 
nearly linear dependence on A. 

Viewed over longer time scales, chaotic ensembles typ- 
ically exhibit a coarse-grained evolution, exponential in 
time, towards a near-uniform population of those easily 
accessible phase space regions not seriously obstructed by 
entropy barriers. For 'wildly chaotic' ensembles, the rate 
A associated with this convergence tends to be insensi- 
tive to the choice of the diagnostic used to probe the 
convergence: different moments and/or coarse-grained 
distribution functions converge at comparable rates. For 
'sticky' ensembles, the convergence can depend more sen- 
sitively on the choice of diagnostic so that, oftentimes, 
one can distinguish between diagnostics characterised by 
'fast' and 'slow' time scales. This behaviour seems in- 
sensitive to the form of the potential. Different 'generic' 
potentials exhibit the same qualitative behaviour; and 
similar behaviour is also observed for 'special' potentials 
like (2.4) where chaotic orbits behave in a near-regular 
fashion most of the time. 

If entropy barriers like cantori or an Arnold web are 
largely irrelevant, this near-invariant distribution may 
correspond (at least approximately) to a true equilib- 
rium, i.e., a uniform, microcanonical population of the 
accessible phase space regions. If, however, such barriers 
arc important, this distribution may differ significantly 
from a true equilibrium. Only over much longer time 
scales will orbits diffuse through the entropy barriers to 
approach a true invariant distribution. 

Viewed over sufficiently long time scales, regular en- 
sembles also exhibit a coarse-grained evolution towards 
an invariant distribution which, now, corresponds to a 
uniform population of the tori to which their motions 
are restricted. This later time evolution is neither strictly 
exponential nor power law, but one can at least identify 
a simply scaling in terms of A, the 'size' of the region 
sampled by the initial conditions. 

Weak perturbations modeled as Gaussian white noise 
can significantly enhance the efficiency of chaotic phase 
mixing, even over time scales sufficiently short that en- 
ergy is almost conserved. In part by relaxing the con- 
straints associated with Liouville's Theorem, the noise 
can help 'fuzz out' short scale structure; and, by facili- 
tating diffusion through entropy barriers, it can dramati- 
cally accelerate the approach towards a true invariant dis- 
tribution. As probed by the rate at which an ensembles 
approach an invariant or near-invariant distribution, the 
efficacy of the noise exhibits a logarithmic dependence on 
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amplitude, the same dependence exhibited by noise as a 
source of accelerated diffusion through entropy barriers. 

Coloured noise with an extremely short autocorrela- 
tion time t c has virtually the same effect as white noise. 
However, when t c becomes comparable to the orbital 
time scale t cr , the effects of the noise begin to decrease; 
and, for t c 3> t cr , coloured noise has only a compara- 
tively minimal effect. For intermediate times, the efficacy 
of coloured noise exhibits a roughly logarithmic depen- 
dence on t c . The detailed form of the noise seems largely 
irrelevant: all that seems to matter are rj and t c . 

The simulations described here also reinforce the in- 
terpretation [jl5| that noise effects chaotic orbits through 
a resonant coupling. Noise can be viewed as a superposi- 
tion of periodic drivings with different frequencies com- 
bined with random phases, so that it is natural to focus 
on how a chaotic ensemble is impacted by driving with a 
single frequency. However, a plot of the efficacy of driv- 
ings with different frequencies, as probed by changes in 
energy induced by the perturbation, shows distinct simi- 
larities to a plot of the power spectra for the unperturbed 
orbits. In particular, the driving has its largest effect for 
frequencies comparable to, or somewhat larger than, the 
frequencies where the spectra have the largest power. 

Given that noise can be viewed as a superposition of 
periodic disturbances with a variety of different frequen- 
cies all combined with random phases, it is not surprising 
that periodic driving impacts orbit ensembles in a fash- 
ion very similar to noise. The only obvious difference 
is that the dependence on frequency is more sensitive 
than the dependence on autocorrelation time exhibited 
by coloured noise. 
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FIG. 1. The dispersion a x for ensembles of regular orbits in 
the two-dimensional dihedral potential with a = 1 and energy 
E = 2.0, generated from 20000 initial conditions sampling a 
configuration space square of side A. The dots correspond 
to data points recorded at intervals St — 4.0. The solid 
curves represent box-car averages over 10 adjacent points. 
The straight line exhibits a best fit approximation to a lin- 
ear growth law. (a) A = 1 x 10~ 3 . (b) A = 5 x 10~ 3 . (c) 
A = 1 x 10~ 2 . (d) A = 2 x 10~ 2 . (e) A = 5 x 10 
best fit growth time to for a linear growth law a 
a function of 'size' A. The dashed line with slope —1 corre- 
sponds to a proportionality to oc 1/A. 
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FIG. 2. Convergence of initially localized ensembles of reg- 
ular orbits towards an invariant distribution. The ensembles 
are the same as those used to generate Fig. 1 d and e, with 
A = 2 x 10~ 2 and A = 5 x 10~ 2 . (a) and (b) The mean value 
(x(t)) and In | <as(*)> | for A = 5 x 10~ 2 . (c) and (d). The same 
for A = 1 x 10~ 2 . (e) and (f) The same as (a) and (b) but 
restricted to a shorter time interval. 




FIG. 3. Smoothed and unsmoothed dispersions a x (a and 
b), (jy (c and d), and o z (e and f) computed for initially 
localised ensembles of 'wildly chaotic' (left panels) and 'sticky' 
chaotic (right panels) orbits with the same energy evolved in 
the ellipsoid plus black hole potential (2.4) 
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FIG. 4. Smoothed dispersions computed for initially lo- 
calised ensembles of chaotic orbits evolved in the potential 
(2.2) with a = b = c = 1. (a) a x (thick solid), a y (solid), and 
<7 Z (dashed) for a 'wildly chaotic' ensemble with E = 1.0. (b) 
cr V x, o~v y , and a vz for the same ensemble, (c) a x , a y , and a z for 
a 'sticky' chaotic ensemble with E — 4.0. (d) a x , a y , and a z 
for a chaotic ensemble with E = 4.0 which is 'wildly chaotic' 
in the x- and y-directions but sticky in the z-direction. 



FIG. 5. Convergence towards a near-invariant distribution 
for an initially localised ensemble of 'wildly chaotic' orbits 
evolved in the potential (2.2) with E — 2.0. (a) Config- 
uration space moments D In (x) (thick solid) , D In (y) + 3 
(thin solid), and D\n(z) + 6 (dotted), (b) Velocity mo- 
ments D\n(v x ) (thick solid), Dln(v x ) +3 (thin solid), and 
D\n(v x ) + 6 (dotted), (c) Mixed moments D\n(xy) (thick 
solid), Dln(yz}+3 (thin solid), and D In (xz)+6 (dotted), (d) 
Dln{xv x ) (thick solid) , D\n(yv y ) (thin solid), and D In 
(dotted), (e) lnD 1 f(x,y) (thick solid), \nD 1 f(v x ,v y ) (thin 
solid), and InD f(x,v x ) (dotted), (f) In D 1 f(x, y) computed 
for different k x k binnings with (from top to bottom) k = 40, 
30, 20, and 10. 
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FIG. 6. The same as the preceding Figure, except now gen- 
erated for a 'sticky' ensemble with E = 4.0. 
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FIG. 7. (a) Representative probes of the rate A of 
convergence towards a near-invariant distribution for a 
'wildly chaotic' ensemble with E = 2.0 evolved in the 
two-dimensional dihedral potential with a = 1. (b) The same 
for a sticky ensemble with E = 4.0. (c) The ratio A/x for the 
rates A associated with D\n{q a ) (diamonds) and In Df (tri- 
angles) for representative 'wildly chaotic ensembles' evolved 
in the dihedral potential with different energies. 
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FIG. 8. (a) A collection of 20000 orbits sampling a 
near-invariant distribution with energy E — 6.0 in the 
two-dimensional dihedral potential with a = 1 at time 
t — 36.0. (b) The same orbits sampled at a much later time 
t = 720. Both ensembles were generated from a localised col- 
lection of initial conditions corresponding to 'wildly chaotic' 
orbits, (c) The same orbits, now evolved in white noise with 
6 = 6.0 and rj = 10~ 6 , again sampled at t = 36.0. (d). White 
noise with O = 6.0 and r\ = 10~ 5 , sampled at t = 36.0. 
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FIG. 9. The quantity ln|(x)| computed for an initially lo- 
calised ensemble of regular orbits with E — 2.0 evolved in 
the potential (2.2) with a = b = c = 1. (a) No noise, (b) 
r) = 10~ 7 . (c) 77 = 10~ 6 . (d) Tj = 10~ 5 . 
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FIG. 10. The dispersion a x computed for initially localised 
ensembles of evolved in the two-dimensional dihedral poten- 
tial with a = 1 in the presence of friction and white noise 
with = E and (from bottom to top) variable r\ = 0, 1CF 7 , 
1(T 5 , and 1CT 4 . The three upper curves were shifted upwards 
by amounts 2, 4, and 6. (a) A 'wildly chaotic' ensemble with 
E = 2.0. (b) A 'sticky' ensemble with E = 4.0. 
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FIG. 11. The quantity Df(x,y) computed for the 'wildly 
chaotic' (a) and 'sticky' (b) ensembles used to generate the 
preceding Figure, again staggered by translations of 2, 4, and 
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FIG. 12. A, the rate at which f(x, y) approaches a 
near-invariant distribution, computed as a function of r) for 
the 'wildly chaotic' ensemble used to generate Fig. 11. The 
dashed line corresponds to the zero noise limit. 
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FIG. 13. (a) The dispersion a x computed for the ensemble 
used to generate Figs. 11 b and 12 b, now integrated in the 
presence of Ornstein-Uhlenbeck noise with = 4.0, i) — 10~ 5 
and (from top to bottom) a = 100.0, a — 10.0, a — 3.33, 
a — 1.0, a = 0.33, and a = 0.1, corresponding to autocorre- 
lation times U = 0.01, 0.1, 0.33, 1.0, 3.33, and 10.0. (b) The 
quantity In computed for the same ensembles. 
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FIG. 14. The analogue of Fig. 13, now computed for orbits 
in the presence of coloured noise characterised by the more 
complex autocorrelation function (2.12). 
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FIG. 15. (a) A, the rate at which f(x, y) approaches a 
near-invariant distribution, computed as a function of t c for 
a 'wildly chaotic' ensemble with E = 2.0 evolved in the pres- 
ence of Ornstein-Uhlenbeck noise with = 1.0 and r\ = 10 -5 . 
The lower dashed line corresponds to evolution in the absence 
of noise. The upper line corresponds to the t c — white noise 
limit, (b) The analogue of (a), now computed for noise with 
the autocorrelation function (2.12). 
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FIG. 16. (a) The maximum value of ge for the ensemble 
used to generate Figs. 10 and lib, now evolved for a time 
t = 256 in the presence of periodic driving of the form given 
by eq. (2.7) with fixed amplitude A = 10~ 2 but different fre- 
quencies u), expressed in units of the original energy Eo = 4.0. 
(b) The same for driving of the form given by eq. (2.5), again 
with A = 10 2 . (c) The composite power spectra, |x(w)| (solid 
curve) and |y(w)| (dashed) computed for the unperturbed en- 
semble. 
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FIG. 17. (a) The dispersion a x computed for an ensemble 
evolved in the two-dimensional dihedral potential with a = 1 
in the absence of all perturbations (lower curve) and in the 
presence of periodic driving of the form (2.7) with A = 10 -1 ' 6 
and (from bottom to top) u> = 0.1, u> = 1.0, u> = 3.0, u> = 10.0, 
and u> = 30.0. The upper curves were shifted upwards by 
amounts 4, 8, 12, 16, and 20. (b) The quantity In for the 
same ensembles, the upper curves again staggered. 



